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SECTION  1 
INTRODUCTION 


Chloral  hydrate  is  one  of  the  oldest  of  synthetic  sedative/hypnotic  drugs  used  in  human 
medicine  (Leibreich,  1869).  Its  popularity  has  decreased  considerably  after  introduction  of 
barbiturates  because  barbiturates  are  much  more  convenient  to  administer.  Chloral  hydrate  must 
be  given  as  a  1  -  2  g  oral  dose  in  a  flavored  solution  to  disguise  its  unpleasant  taste,  or  as  a 
retention  enema  or  suppository.  Despite  administration  in  huge  doses,  its  toxicity  is  very  low 
(lethal  dose  of  about  10  g  per  70  kg  man)  and,  typically,  its  only  side  effect  is  some  gastric 
irritation  after  therapeutic  oral  administration  (Goth,  1964).  In  the  organism,  chloral  hydrate  is 
rapidly  metabolized*  mainly  to  trichloroethanol  (TCOH)  and  its  glucuronide  (TCOG)  which  is 
excreted  with  the  urine  (Butler,  1948).  Both  chloral  (CH)  and  TCOH  can  exert  hypnotic  activity 
on  the  central  nervous  system  (MacKay  and  Cooper,  1962). 

About  11%  of  the  administered  chloral  hydrate  is  oxidized  to  trichloroacetic  acid  (TCA) 
(Cabana  and  Gessner,  1970;  Muller  et  al.,  1974),  and  possibly  to  a  small  amount  of 
dichloroacetic  acid  (DCA)  (Davidson  and  Beliles,  1991).  However,  it  is  not  clear  if  DCA  is  a 
metabolic  product  of  CH  as  suggested  by  Dekant  et  al.  (1984),  a  product  of  metabolic  reductive 
dechlorination  of  TCA  as  suggested  by  Larson  and  Bull  (1992),  or  rather  it  is  an  analytical 
artifact  produced  by  interconversion  of  TCA  in  biological  samples  (Ketcha  et  al.,  1995). 

A  realization  that  CH  is  generated  as  a  chlorination  by-product  in  several  municipal 
drinking  water  supplies,  as  well  as  that  it  is  produced  in  the  organism  as  an  intermediate 
metabolite  from  the  ubiquitous  environmental  pollutant  -  trichloroethylene  (TCE)  (Cole  et  al., 
1975;  Reynolds  and  Moslen,  1981)  rejuvenated  interest  in  the  toxicology  of  CH  and  its  hydrate 
(Waters  et  al.,  1977).  Especially,  the  recent  reports  about  mutagenic  and  hepatocarcinogenic 
potential  of  chloral  hydrate  (Daniel  et  al.,  1992)  have  put  its  importance  in  a  new  light.  Chloral 
and  TCE  occur  as  environmental  contaminants  at  Air  Force,  Navy,  and  Army  installations  and 
for  some  time  these  chlorinated  compounds  have  been  part  of  the  TriService  research  initiative. 
For  a  rational  characterization  of  the  potentially  harmful  effects  of  CH  and  its  metabolites,  it  is 
crucial  to  be  able  to  reconstruct  its  distribution  and  conversion  to  its  main  metabolites:  TCOH, 
TCOG,  TCA,  and  DCA,  from  the  CH  exposure  dosage.  The  approach  ideally  suited  for  this 
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purpose  is  the  use  of  a  physiologically  based  pharmacokinetic  (PBPK)  model  (Yang  and 
Andersen,  1994).  This  report  describes  development  of  a  PBPK  model  of  CH  and  its  main 
metabolites,  mathematical  description  of  the  interlinking  between  sub-models  for  each 
metabolite,  model  verification  with  data  from  the  literature  and  its  experimental  calibration  in 
B6C3F1  mice. 
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SECTION  2 

MATERIALS  AND  METHODS 

A  PBPK  model  was  written  in  ACSL,  a  Fortran-based  continuous  simulation  language 
(Mitchell  and  Gauthier,  1987),  and  simulations  were  performed  using  a  SIMUSOLV  software 
package  with  optimization  capabilities  (Steiner  et  al.,  1990)  on  a  VAX/VMS  mainframe 
computer  (VAX8530,  Digital  Equipment  Corp.,  Maynard,  MA).  Parameters  were  optimized  by 
SIMUSOLV  which  uses  the  log  likelihood  function  as  the  criterion,  and  either  the  generalized 
reduced  gradient  method  for  single  parameter  optimization  or  the  Nelder-Mead  search  method 
for  multiple  parameters  optimization  to  adjust  the  values. 

A  method  for  nonvolatile  chemicals,  after  Jepson  et  al.  (1994),  was  used  for  measuring 
tissue  partition  coefficients  for  CH  and  TCOH.  Briefly,  0.5  g  of  blood  or  other  tissue  was  added 
to  a  5  mL  of  solution  containing  CH  or  TCOH  in  an  appropriate  concentration,  20%  of  lead 
acetate,  and  0.9%  NaCl  in  20  mL  vial  capped  with  teflon/rubber  septum.  The  vials  were 
equilibrated  for  18  hr  at  37  °C  with  vortexing  at  a  medium  speed.  Equilibrated  tissue 
supernatants  were  centrifuged  at  1500  rpm  for  10  min  and  the  resultant  supernatants  were 
filtered  through  prewashed  Milipore  filters  (Ultra-PF,  low-binding  cellulose,  10,000  NMWL). 
The  ultrafiltrate  was  extracted  with  ethyl  acetate  and  analyzed  for  CH  or  TCOH  by  gas 
chromatography. 

B6C3F1  male  mice  (30-40  g  body  weight;  Charles  River,  Inc.)  were  dosed  i.v.  (by  the  tail 
vein)  with  chloral  hydrate  or  TCOH  at  an  appropriate  concentration  (10  or  100  mg/kg), 
dissolved  in  physiological  saline.  The  injection  interval  for  each  mouse  was  timed.  At 
appropriate  times,  the  mice  were  sacrificed  with  C02,  and  samples  of  blood  and  liver  were 
collected  to  pre- weighed  vials  containing  20%  lead  acetate.  The  samples  were  homogenized, 
extracted  with  ethyl  acetate,  and  analyzed  for  CH  or  TCOH  by  gas  chromatography. 

The  animals  used  in  this  study  were  handled  in  accordance  with  the  principles  stated  in  the 
"Guide  for  the  Care  and  Use  of  Laboratory  Animals"  prepared  by  the  Committee  on  Care  and 
Use  of  Laboratory  Animal  Resources,  National  Research  Council,  Department  of  Health  and 
Human  Services,  National  Institutes  of  Health,  Publication  No.  86-23,  1985;  and  the  Animal 
Welfare  Act  of  1966,  as  amended. 
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The  tissue  samples  were  analyzed  for  CH  and/or  TCOH  by  a  Hewlett  Packard  5890  Series 
II  gas  chromatograph  equipped  with  a  7673  A  autosampler.  The  blank  and  reference  vials 
contained  no  tissue  but  were  processed  in  the  same  manner  as  the  vials  that  contained  actual 
tissue  sample.  The  gas  chromatograph  conditions  that  separated  CH  from  TCOH  were:  injection 
temp.  =  175  °C,  Electron  Capture  Detector  temp.  =  300  °C;  oven  program  (80  °C  for  4  min; 
25  °C/min  to  180  °C  for  2.2  min);  Vocol  30  M  x  0.53  mm  column.  Data  were  collected  using 
the  Turbochrom  Nelson  Analytical  System  v.4.03.  Standard  curves  were  made  in  the 
appropriate  biological  matrix,  and  treated  with  the  same  procedures  as  test  samples. 

Chloral  hydrate  and  TCOH  were  obtained  from  Sigma  Chemical  Co.,  Ltd.  and  lead  acetate 
used  to  inhibit  metabolism  in  analytical  samples  was  obtained  from  Mallinckrodt.  All  other 
commercial  reagents  used  were  of  analytical  purity. 
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SECTION  3 
RESULTS 

PBPK  Model  Structure 

Figure  1  shows  a  general  scheme  of  an  interlinked  PBPK  model  for  CH  and  its  metabolites. 
The  main  model  consisted  of  three  sub-models  for  CH,  TCOH,  and  TCOG  (treated  as  separate 
objects)  linked  by  rates  of  production  of  TCOH  from  CH  (RAMCH  [mg/hr]),  glucuronide 
production  from  TCOH  (RGLUC  [mg/hr]),  and  reabsorption  of  TCOH  from  hydrolyzed  TCOG 
(RAGLUC  [mg/hr]).  Additional  outputs  from  these  sub-models  were:  the  amount  of  DC  A 
produced  from  CH  (ADC A  [mg])  with  the  rate  RIDA  [mg/hr],  the  amount  of  TCA  (ATCA 
[mg])  produced  from  CH  with  the  rate  RTA  [mg/hr]  and  from  TCOH  with  the  rate  RAMOCH 
[mg/hr],  and  the  amount  of  TCOH  glucuronide  (AGU  [mg])  excreted  to  the  urine  with  the  rate 
RGU  [mg/hr] .  The  output  from  each  sub-model  was  available  for  input  in  subsequent  operations 
of  this  object-oriented  program. 

Rates  of  Mass  Transfer  Between  Sub-models 

i.  Production  of  TCOH  from  CH: 

RAMCH  =  (VMCHOH*CVCHL)/(KMCHOH +CVCHL)  +  PCTCOH*CVCHL*VL 
where:  VMCHOH  [mg/hr]  is  a  pseudo-maximum  velocity  of  TCOH  formation  from  CH, 
adjusted  allometrically  to  the  0.7  power  of  the  body  weight  [kg];  *  is  a  multiplication;  CVCHL 
[mg/L]  is  a  CH  venous  blood  concentration  leaving  the  liver;  KMCHOH  [mg/L]  is  an  apparent 
Michaelis  constant  of  TCOH  formation  from  CH;  PCTCOH  [1 /hr/animal]  is  a  first  order  rate 
constant  of  TCOH  formation  from  CH,  adjusted  allometrically  to  the  0.3  power  of  the  body 
weight  [kg];  VL  [L]  is  a  volume  of  liver,  adjusted  allometrically  to  the  body  weight. 

ii.  Production  of  TCOG  from  TCOH: 

RGLUC  =  (VMTCOH*CVOHL)/(KMTCOH+CVOHL) 
where:  VMTCOH  [mg/hr]  is  a  pseudo-maximum  velocity  of  TCOG  formation  from  TCOH, 
adjusted  allometrically  to  the  0.7  power  of  the  body  weight  [kg];  CVOHL  [mg/L]  is  a  TCOH 
venous  blood  concentration  leaving  the  liver;  KMTCOH  [mg/L]  is  an  apparent  Michaelis 
constant  of  TCOG  formation  from  TCOH. 
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Interlinked  PBPK  model  for  CH  and  its  metabolites 
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iii.  Reabsorption  of  TCOH  from  hydrolyzed  glucuronide: 

RAGLUC  =  KAOHBI*ATGB*eKAOHBI*T 

where:  KAOHBI  [1/hr]  is  a  TCOH  uptake  from  gut  rate  constant;  ATGB  [mg/animal]  is  the 
amount  of  TCOH  formed  by  complete  hydrolysis  of  glucuronide  in  the  intestine;  T  [hr]  is  the 
time. 

Rates  of  Metabolite  Mass  Output  from  Sub-models 

i.  Formation  of  DCA  from  CH: 

RIDA  =  CVCHL*PIDA*VL 

where:  CVCHL  [mg/L]  is  a  CH  venous  blood  concentration  leaving  the  liver;  PIDA 
[1/hr/animal]  is  a  first  order  rate  constant  of  DCA  formation  from  CH,  adjusted  allometrically 
to  the  0.3  power  of  the  body  weight  [kg];  VL  [L]  is  a  volume  of  liver,  adjusted  allometrically 
to  the  body  weight. 

ii.  Formation  of  TCA  from  CH: 

RTA  =  CVCHL*PCTCA*VL 

where:  PCTCA  [1  /hr/animal]  is  a  first  order  rate  constant  of  TCA  formation  from  CH,  adjusted 
allometrically  to  the  0.3  power  of  the  body  weight  [kg]. 

iii.  Formation  of  TCA  from  TCOH: 

RAMOH  =  CVOHL*PCOTA*VL 

where:  PCOTA  [1/hr/animal]  is  a  first  order  rate  constant  of  TCA  formation  from  TCOH, 
adjusted  allometrically  to  the  0.3  power  of  the  body  weight  [kg]. 

At  the  input  to  the  next  sub-model,  the  rates  (RIDA,  RTA,  and  RAMOH)  were  adjusted 
for  molecular  weight  of  metabolites. 

iv.  Excretion  of  TCOG  with  the  urine: 

RGU  =  CBVG*CLGLU 

where:  CBVG  [mg/L]  is  a  concentration  of  glucuronide  in  ultrafiltrated  blood;  CLGLU  [L/hr] 
is  a  renal  clearence  of  TCOG  adjusted  allometrically  to  the  body  weight  [kg] . 
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Compartment  Description  and  Governing  Equations 

The  basic  assumption  in  constructing  this  PBPK  model  was  that  blood  flow  to  the  tissue  is 
limiting  the  metabolite  "m"  delivery.  Because  any  metabolite  is  retained  by  the  tissue  according 
to  its  tissue/blood  partition  coefficient  (PJ)  which  may  be  measured  in  vitro,  the  concentration 
of  the  metabolite  in  venous  blood  leaving  the  tissue  (CVTO)  during  the  equilibration  phase  is 
lower  than  the  concentration  in  arterial  blood  (CmA).  Therefore,  the  rate  of  change  of  metabolite 
amount  in  tissue  (dA^/dt),  is  given  by  a  simple  difference  between  concentration  in  blood 
entering  and  exiting  the  tissue  (C^A-CV^)  multiplied  by  the  blood  flow  (Q,). 

Integrating  this  equation  over  a  given  time,  one  can  calculate  the  amount  of  metabolite 
present  in  tissue  (Ami)  and,  therefore,  if  the  actual  volume  of  tissue  (V()  is  known,  one  can 
calculate  the  concentration  of  substance  in  the  tissue  (Cmi)  at  any  time.  Using  these  simple 
principles,  the  mass  transfer  equations  describing  each  compartment  building  the  sub-models  for 
CH,  TCOH,  and  TCOG  (schematically  shown  in  Figures  2  and  3)  are  defined  below. 

For  each  well-stirred  compartment  without  metabolism  or  other  losses  (fat  tissue,  gut, 
slowly  perfused  and  rapidly  perfused  tissues),  the  rate  of  change  in  the  amount  "A"  [mg]  of 
metabolite  "m"  (dAmi  over  time  [hr])  was  defined  as  follows: 

dAm;/dt  =  Qi(CmA  -  CVm() 

where:  subscript  i  represents  "i-th"  compartment;  m  represents  the  metabolite  (CH,  TCOH,  or 
TCOG);  Q,  [L/hr]  represents  the  blood  flow  through  the  i-th  compartment;  C^A  [mg/L] 
represents  the  arterial  concentration  of  metabolite  m;  CVmi  [mg/L]  represents  the  venous 
concentration  of  metabolite  m  leaving  the  i-th  compartment  (CVTO  =  Cm(/PTO;  where  Cmi  [mg/kg] 
is  a  concentration  of  metabolite  m  in  the  tissue  in  i-th  compartment  and  Pm,  is  the  tissue/blood 
partition  coefficient  of  metabolite  m  for  i-th  compartment.  Cmi  =  ATO/V„  where  V,  [kg] 
represents  the  volume  of  the  i-th  compartment). 

For  the  liver  compartment  a  loss  term  (RAMm  [mg/hr])  was  added  to  the  well-stirred 
compartment  description  to  account  for  metabolism,  and  an  increment  term  RAMmI  [mg/hr]  to 
account  for  production  in  the  previous  sub-model  (RAMm  and  RAMmI  are  the  rates  of  mass 
transfer  between  sub-models  or  the  rates  of  mass  output  from  sub-models,  adjusted  for  molecular 
weight  of  the  substrate).  As  explained  above,  depending  on  the  metabolite,  rates  of  mass  transfer 
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or  mass  output  followed  either  the  Michaelis-Menten  kinetic  equation,  mixed  Michaelis-Menten 
and  first  order  rate  kinetic  equation,  or  the  first  order  rate  of  metabolism: 

dA^L/dt  =  QL(CmA  -  CVmL)  -  RAMm  +  RAMmI 

Because  both  CH  and  TCOH  (but  not  TCOG)  are  slightly  volatile,  for  the  lung 
compartment  with  two  theoretically  possible  mass  inputs  (mixed  venous  blood  and  inhaled  air) 
and  two  theoretically  possible  outputs  (arterial  blood  and  exhaled  air),  at  steady  state  the  amount 
in  alveolar  air  is  in  equilibrium  with  the  amount  in  lung  blood.  Therefore: 

QP(CIm  -  CXJ  =  QC(CmA  -  CmV) 

CXm  =  CmA/PmB 

where:  QP  [L/hr]  is  the  air  flow  through  the  lungs  (alveolar  ventilation  rate),  CIm  [mg/L]  is  the 
concentration  in  inhaled  air,  C^X  [mg/L]  is  the  concentration  in  alveolar  air,  C„A  [mg/L]  is  the 
arterial  concentration  (leaving  the  lungs),  PmB  is  blood/air  partition  coefficient,  QC  [L/hr]  is  the 
blood  flow  through  the  lungs  (rate  of  cardiac  output),  CmV  [mg/L]  is  the  venous  concentration 
(entering  the  lungs).  This  equation  is  solved  for  C^A. 

The  rate  of  change  in  the  amount  of  CH  or  TCOH  (RMRm  [mg/hr])  in  the  gastrointestinal 
tract  (after  single  gavage  dosing)  was  defined  analogously  to  the  reabsorption  of  TCOH  from 
hydrolyzed  glucuronide  in  the  gut  (as  described  above): 

RMRm  =  -KAm*DOSm*e-KAm*T 

where:  KAm  [l/hr]  is  a  CH  or  TCOH  uptake  from  gut  rate  constant;  DOSm  [mg/animal]  is  the 
total  oral  dose;  T  [hr]  is  the  time. 

The  rate  of  change  in  the  amount  of  TCOG  remaining  in  the  body  (RBODYG  [mg/hr])  was 
defined  as  a  difference  between  the  rate  of  retention  or  wash-out  from  the  tissues  and  the  amount 
excreted  with  the  urine: 

RBODY  =  QB*(CGB-CBVG)  -  RGU 

where:  QB  [L/hr]  is  a  blood  flow  through  the  body  except  liver;  CGB  [mg/L]  is  a  concentration 
of  TCOG  in  mixed  blood;  CBVG  [mg/L]  is  a  concentration  of  TCOG  in  ultrafiltrated  blood  (in 
equilibrium  with  the  ultrafiltrate);  RGU  [mg/hr]  is  the  rate  of  excretion  of  TCOG  with  the  urine 
(defined  above). 
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exhaled  air  [mg/L];  COHA  -  concentration  of  TCOH  in  arterial  blood  [mg/L];  CGBIL  -  concentration  of  TCOG 
in  the  bile  [mg/L];  CVOHG  -  concentration  of  TCOH  in  portal  vein  [mg/L];  CBVG  -  concentration  of  TCOG  in 
mixed  venous  blood  [mg/L];  AMOH  -  amount  of  TCOH  metabolized  [mg];  QP  -  alveolar  ventilation  [L/hr]; 
COHV  -  concentration  of  TCOH  in  mixed  venous  blood  [mg/L];  IVDOH  -  intravenous  dose  of  TCOH  [mg/kg]. 


Numerical  Values  of  PBPK  Model  Constants  Used  in  Simulations 

Physiological  parameters  used  in  the  PBPK  model  for  simulations  were  adapted  from  the 
literature  (Compilation  by  Lindstedt,S.:  Unpublished  physiological  parameters.  Physiological 
Parameters  Working  Group,  ILSI  Risk  Science  Institute;  and  Arms  and  Travis,  1988).  These 
parameters  included  alveolar  ventilation  (QPC  [L/hr/kg]),  cardiac  output  (QCC  [L/hr/kg]),  tissue 
blood  flows  (Q.  [fraction  of  QCC]),  tissue  volumes  (V,C  [fraction  of  body  weight]),  and  are 
listed  in  Table  1. 


TABLE  1.  Physiological  Parameters  Used  for  PBPK  Model  Simulations 


Parameter 

Description 

Value 

[Unit] 

BW 

Body  weight 

Measured 

[kg] 

QPC 

Alveolar  ventilation 

30.0 

[L/hr/kg] 

QCC 

Cardiac  output 

16.5 

[L/hr/kg] 

QGC 

Blood  flow  to  gut 

0.175 

[ratio] 

QLC 

Blood  flow  through  hepatic  artery 

0.24  -  QGC 

[ratio] 

QFC 

Blood  flow  to  fat 

0.05 

[ratio] 

QSC 

Blood  flow  to  slowly  perfused  tissues 

0.238 

[ratio] 

QRC 

Blood  flow  to  rapidly  perfused  tissues 

0.472 

[ratio] 

QKC 

Blood  flow  to  kidney 

QRC  -  0.252  [ratio] 

QUC 

Urine  flow 

0.0006 

[L/hr/kg] 

QBILC 

Bile  flow 

0.00015 

[L/hr/kg] 

VLC 

Volume  of  liver 

0.05 

[ratio] 

VFC 

Volume  of  fat 

0.04 

[ratio] 

VSC 

Volume  of  slowly  perfused  tissues 

0.558 

[ratio] 

VRC 

Volume  of  rapidly  perfused  tissues 

0.031 

[ratio] 

VGC 

Volume  of  gut  tissue 

0.033 

[ratio] 

VKC 

Volume  of  kidney 

0.018 

[ratio] 

VBLD 

Volume  of  blood 

0.06 

[ratio] 

Chloral  hydrate  and  TCOH  are  very  soluble  in  blood  so  their  partitioning  between  blood 
and  air  could  not  be  measured  with  the  method  employed.  Therefore,  their  blood/air  partition 
coefficients  were  estimated  from  solubilities  of  radioactive  analogues  of  similar  compounds 
reported  in  the  literature.  Similarly,  partitioning  of  TCOG  in  bile,  blood,  and  solid  tissues,  as 
well  as  renal  clearance,  were  estimated  from  the  glucuronide  distribution  data  available  in  the 
literature  (Garrett  and  Lambert,  1973).  Partition  coefficients  of  CH  and  TCOH  in  other  tissues 
were  measured  in  our  laboratory  (Seckel  et  al.,  1995).  These  physicochemical  parameters  are 
listed  in  Table  2.  Metabolic  parameters  and  constants  for  conversions  of  CH  to  TCOH,  TCOG, 
TCA,  and  DCA  were  initially  estimated  from  the  literature  (Cabana  and  Gessner,  1970;  Garrett 
and  Lambert,  1973;  Larson  and  Bull,  1992;  Templin  et  al.,  1993)  and  later,  fitted  and  optimized 
with  the  PBPK  model.  The  final  set  of  metabolic  parameters  is  listed  in  Table  3. 
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TABLE  2.  Physicochemical  Parameters  Used  for  PBPK  Model  Simulations 


Parameter 

Description 

Value 

[Unit] 

Partition  coefficients  for  chloral 

hydrate 

PCHB 

Estimated  blood/air 

500.0 

[ratio] 

PCHL 

Measured  liver/blood 

1.47 

[ratio] 

PCHF 

Measured  fat/blood 

0.48 

[ratio] 

PCHS 

Measured  slowly  perfused/blood 

1.35 

[ratio] 

PCHR 

Measured  rapidly  perfused/blood 

1.47 

[ratio] 

PCHG 

Measured  gut/blood 

1.47 

[ratio] 

Partition  coefficients  for  trichloroethanol 

POHB 

Estimated  blood/air 

5000.0 

[ratio] 

POHL 

Measured  liver/blood 

1.06 

[ratio] 

POHF 

Measured  fat/blood 

1.53 

[ratio] 

POHS 

Measured  slowly  perfused/blood 

1.11 

[ratio] 

POHR 

Measured  rapidly  perfused/blood 

1.06 

[ratio] 

POHG 

Measured  gut /blood 

1.06 

[ratio] 

Partition  coefficients  for  trichloroethanol  glucuronide 

PCBIL 

Estimated  bile/blood 

70.0 

[ratio] 

PCGBO 

Estimated  body/blood 

0.31 

[ratio] 

CLGLUC 

Estimated  renal  clarence 

0.65 

[L/hr/kg] 

Molecular  weights 

MWCH 

Molecular  weight  CH 

147.5 

[g/mol] 

MWTCOH 

Molecular  weight  TCOH 

149.5 

[g/mol] 

MWGLUC 

Molecular  weight  TCOG 

325.4 

[g/mol] 

MWTCA 

Molecular  weight  TCA 

163.5 

[g/mol] 

MWDCA 

Molecular  weight  DCA 

129.0 

[g/mol] 

TABLE  3. 

Metabolic  Parameters  Used  for  PBPK  Model  Simulations 

Parameter 

Description 

Value 

[Unit] 

PCCOH 

First  order  metabolic  constants 

Conversion  CH->TCOH 

15.0 

[l/hr/kg] 

PCCA 

Conversion  CH->TCA 

0.02 

[l/hr/kg] 

PCIDA 

Conversion  CH->DCA 

0.3 

[l/hr/kg] 

KCHC 

Metabolic  loss  of  CH  by  blood 

1.04 

[l/hr/kg] 

PCOA 

Conversion  TCOH->TCA 

3.5 

[l/hr/kg] 

KOHDC 

Metabolic  loss  of  TCOH  by  blood 

7.9 

[l/hr/kg] 

VMCHOC 

Michaelis-Menten  metabolic  constants 
Pseudo-V^  CH->TCOH 

3.23 

[mg/hr/kg] 

VMTCOC 

Pseudo-V^  TCOH-  >TCOG 

15.0 

[mg/hr/kg] 

KMCHOH 

Apparent  Michaelis  constant  CH->TCOH 

0.0675 

[mg/L] 

KMTCOH 

Apparent  Michaelis  constant  TCOH->TCOG 

1.0 

[mg/L] 

Experimental  Calibration  of  PBPK  Model 

Before  the  experimental  calibration  in  mice,  the  PBPK  model  was  verified  with  available 
literature  data.  The  results  of  PBPK  model  computer  simulations  of  historical  experimental  data 
by  Cabana  and  Gessner  (1970)  from  mice  treated  i.v.  with  500  mg  of  CH  per  kg,  are  shown  in 
Figure  4.  The  data  set  for  body  levels  of  remaining  CH,  shown  in  Figure  4a,  was  used  to  adjust 
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initial  PBPK  model  parameters.  Once  the  parameters  affecting  directly  rates  of  mass  output  from 
the  CH  sub-model  (to  TCOH,  TCA,  and  DCA)  were  optimized  to  fit  available  literature  data 
(Figure  4a;  Larson  and  Bull,  1992;  Templin  et  al.,  1993),  the  other  data  from  Cabana  and 
Gessner  (1970)  for  body  levels  of  total  TCOH  (Figure  4b),  and  free  (non-glucuronidated)  TCOH 
(Figure  4c)  were  poorly  predicted  by  the  PBPK  model  with  interlinking  rate  constants  (PCCOH, 
VMCHOC  and  KMCHOH)  optimized  to  fit  experimental  points  from  Figure  4a.  Therefore,  the 
coupling  metabolic  rate  constants  linking  the  CH  and  TCOH  sub-models  (affecting  directly 
RAMCH)  were  estimated  as  a  compromise  between  values  giving  statistically  "the  best"  fit  to 
experimental  points  from  Figure  4a,  b,  and  c.  These  values  are  listed  in  Table  3,  and  the  results 
of  computer  simulation  are  plotted  on  a  semi-logarithmic  scale  as  continuous  lines  in  Figures  4a, 
b,  and  c. 

It  was  not  possible  to  find  in  the  literature  a  coherent  experimental  data  set  describing 
TCOH  and  TCOG  distribution  in  mice,  suitable  for  PBPK  simulation.  Because  of  experimental 
constraints,  it  was  especially  difficult  to  find  experimental  results  on  excretion  and  reabsorption 
from  the  bile  in  mice  or  other  small  rodents.  However,  Garrett  and  Lambert  (1973)  have 
published  very  detailed  experimental  data  on  distribution  of  TCOH  and  TCOG  in  dogs.  The 
physiological  parameters  used  for  the  dog  were  scaled  allometrically  from  those  presented  in 
Table  1  for  the  mice.  Figures  5  and  6  show  the  results  of  PBPK  model  computer  simulations  of 
these  experimental  data  by  Garrett  and  Lambert  (1973)  from  dogs  treated  i.v.  with  TCOH  or 
TCOG.  The  data  set  for  TCOH  blood  concentrations  in  dogs  treated  with  30  mg/kg  of  TCOH, 
shown  in  Figure  5a,  was  used  to  adjust  the  initial  model  parameters.  Thus,  we  have  checked 
how  well  the  PBPK  model  simulates  distribution  of  a  very  high  dose  of  TCOH  (Figure  5b),  and 
finally,  we  have  adjusted  initial  model  parameters  for  TCOG  (Figure  6a).  Once  the  parameters 
directly  affecting  rates  of  mass  output  from  the  TCOH  sub-model  (to  TCA  and  TCOG)  were 
optimized  to  fit  the  available  literature  data  (Figure  5a),  the  other  data  set  from  Garrett  and 
Lambert  (1973)  for  simultaneous  blood  concentrations  of  TCOH  (COHV  [mg/L]),  TCOG  (CGB 
[mg/L]),  and  the  concentration  of  TCOG  in  the  bile  (CGBIL  [mg/L])  was  simulated  without 
changing  the  previously  estimated  parameters  (Figure  6b).  These  parameters  are  listed  in 
Table  3. 

Using  the  complete  set  of  physiological,  physicochemical,  and  metabolic  parameters 
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calibrated  with  data  from  the  literature  (Tables  1,  2,  and  3),  we  have  simulated  our  own 
experimental  data  from  B6C3F1  mice  treated  i.v.  with  10  or  100  mg/kg  of  chloral  hydrate  or 
TCOH  (Figure  7).  Each  experimental  time-point  corresponds  to  the  sample  collected  from  one 
animal. 


Pharmacokinetics  of  Trichloroethanol  in  Dogs 

TCOH  iv  30n,g/k9  00G  (Garrett, 1973)  TCOH  iv  BOmg/kg  DOGIGarrett, 1973) 
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Pharmacokinetics  of  TCOH  Glucuronide  in  Dogs 

TCOG  iv  450mg/kg  00G( Garrett,  1973)  TCOH  iv  1 1  rug/Kg  DOGIGarret-t,  1973) 
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Pharmacokinetics  of  Chloral  and  Trichloroethanol  in  Mice 
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SECTION  4 
DISCUSSION 

The  developed  interlinked  PBPK  model  for  CH  and  its  metabolites  represents  an  attempt 
to  apply  an  "object  oriented"  programming  strategy  (Moniz  Barreto  et  al.,  1994)  in  PBPK 
modeling.  The  main  model  consisted  of  three  sub-models  for  CH,  TCOH,  and  TCOG.  Each  sub¬ 
model  was  calibrated  separately  (as  the  separate  programming  "object"),  validated  with 
experimental  data,  and  may  be  used  to  simulate  distribution  of  each  metabolite  (Figures  4a,  5a, 
5b,  6a,  7a,  and  7b).  When  linked  by  rates  of  mass  transfer  between  sub-models,  the  interlinked 
PBPK  model  was  validated  for  sequential  production  and  distribution  of  the  metabolic  products 
resultant  from  the  previous  metabolite,  and  may  be  used  to  simulate  distribution  of  each 
metabolite  derived  from  the  parent  compound  (Figures  4b,  4c,  and  6b).  This  approach  may  be 
especially  useful  for  pharmacokinetic  exposure  characterization  and  risk  assessment  when 
metabolic  products  are  more  toxic  than  their  parent  compound.  This  seems  to  be  the  case  with 
TCE  and  CH  (Davidson  and  Beliles,  1991).  The  existing  PBPK  models  for  TCE,  described  in 
the  literature  (Andersen  et  al.,  1987;  Dallas  et  al.,  1991),  addressed  only  the  distribution  of  the 
parent  compound  (Buben  and  O’Flaherty,  1985).  The  PBPK  model  developed  by  Fisher  et  al. 
(1989;  1990;  1991),  however,  described  also  the  production  of  TCE  metabolite,  TCA.  Our 
model  was  based  on  more  detailed  mechanistic  information  about  production  and  disposition  of 
intermediate  metabolites,  CH,  TCOH,  and  TCOG,  which  have  distinct  pharmacodynamic 
properties,  different  than  TCA  (Davidson  and  Beliles,  1991).  Thus,  our  attempt  to  interlink  the 
sub-models  into  one  functional  PBPK  model  describing  the  production  of  several  metabolites, 
may  be  considered  as  a  novel  but  consequent  endeavor  in  a  pharmacokinetic  description  of  TCE 
and  its  main  metabolites. 

Each  of  our  sub-models  was  constructed  in  accordance  with  the  conventional  flow  rate 
limited  PBPK  modeling  routine  (Yang  and  Andersen,  1994)  and  could  be  allometrically  scaled 
according  to  the  body  weight  of  the  animal.  As  demonstrated  in  Figures  4b,  5a,  and  7b,  this 
allowed  us  to  describe  successfully  the  pharmacokinetics  of  TCOH  in  as  distant  animal  species 
as  a  mouse  and  a  dog.  Obviously,  the  pharmacokinetic  model  may  be  only  as  good  as  are  the 
data  used  for  its  calibration  and  validation.  Since  we  were  unable  to  identify  in  the  available 
literature  a  coherent  data  set  for  pharmacokinetics  of  chloral  hydrate  and  its  main  metabolites 
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in  humans,  it  is  still  to  be  proved  that  this  model  can  also  adequately  describe  the  distribution 
of  CH  in  men.  On  the  other  hand,  PBPK  models  based  on  the  same  mathematical  and 
physiological  principles  were  already  used  successfully  to  estimate  the  exposure  to  other 
chemicals  in  humans  under  controlled  (e.g.,  Byczkowski  and  Fisher,  1994)  and  environmental 
scenarios  (e.g.,  Fisher  and  Allen,  1993;  Byczkowski  and  Fisher,  1995).  Therefore,  it  seems  that 
once  adequately  validated  in  humans,  the  developed  PBPK  model  might  be  applied  not  only  for 
a  computer-aided  simulation  of  body  levels  of  chloral  hydrate  in  a  therapeutic  situation,  but  also 
for  the  estimate  of  toxicokinetics  of  its  active  metabolites  generated  during  the  environmental 
pollution  scenario. 
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APPENDIX 

Source  Codes  of  BBPD  Model  Written  in  ACSL.  The  *.CSL  and  *.CMD  files  should  be 
executed  under  SIMUSOLV. 

*.CSL  File 


PROGRAM:  INTERLINKED  PBPK  MODEL  -  CHLORAL  AND  ITS  MAIN  METABOLITES  - 
'NOTES:  codes  by  Janusz  Z.  Byczkowski  03/04/94 , revs . 09/12/94 , 02/07/95 


'1)  ooPBPK  CH  AND  ITS  METABOLISM  SUB-MODEL  -  DATA  FROM  MOUSE, 

'  Janusz  Z.  Byczkowski  02/14/94,  rev.  02/07/95,  verified  with 

'  historical  data:  mouse  i.v.-  Cabana  &  Gessner  (1970)  06/28/94, 

'  calibrated:  mouse  i.v.-  Seckel  et  al.  10/20/94, 

'  PCs  from  mouse:  08/30/94,  and  Seckel  et  al.  10/20/94 
'2)  ooPBPK  TCOH  AND  ITS  METABOLISM  SUB -MODEL  -  DATA  FROM  DOG 
'  AND  MOUSE, 

'  Janusz  Z.  Byczkowski  02/16/94,  rev.  02/07/95,  verified  with 
'  historical  data:  dog  i.v.-  Garrett  &  Lambert  (1973),  04/05/94 
'  calibrated:  mouse  i.v.-  Seckel  et  al.  11/17/94, 

'  PCs  from  mouse  08/31/94,  and  Seckel  et  al .  10/20/94 
'3)  ooPBPK  TCOG  SUB -MODEL  -  DATA  FROM  DOG,  AND  MOUSE, 

'  Janusz  Z.  Byczkowski  02/18/94,  rev.  02/07/95,  verified  with 

'  historical  data:  dog  i.v.  -  Garrett  &  Lambert  (1973),  04/20/94 

'  calibrated:  mouse  i.v.-  Seckel  et  al .  11/17/94, 

'  PCs  for  TCOG  estimated  from  dog  04/25/94 

'  OUTPUT  OF  TCA  and  DCA  verified  with  historical  data:  mouse 

'  Larson  &  Bull  (1992)  and  Templin  et  al.  (1993)  09/12/94. 

INITIAL 

'SCALED  TO  MICE 

LOGICAL  CC  $'Flag  set  to  .TRUE,  for  closed  chamber  runs 

'  ***  PHYSIOLOGICAL  PARAMETERS 


'ANIMAL  PARAMETERS 
CONSTANT  BW  =  0.03 

CONSTANT  QPC=  30. 

CONSTANT  QCC=  16.5 

CONSTANT  QGC=  0.175 

QLC=0 . 24-QGC 


CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 


QFC= 

QSC= 

QRC= 

QUC= 


0.05 

0.238 

0.472 


$' 

$' 

$' 

$' 

$' 

$' 

$' 

$' 


0.0006$ 


QBILC=0. 00015$' 


QKC=QRC- 0 .252  $ 


Body  weight  (kg) 

Alveolar  ventilation  rate  (1/hr) 

Cardiac  output  (1/hr) 

Fractional  blood  flow  to  gut 

Fractional  blood  flow  through  hepatic  artery 

Fractional  blood  flow  to  fat 

Fractional  blood  flow  to  slow 

Fractional  blood  flow  to  rapid 

Urine  flow  (L/hr*kg-l) 

Bile  flow  (L/hr*kg-l) 

Fractional  blood  flow  to  kidneys 


CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 

CONSTANT 


VLC=  0.05 
VFC=  0.10 
VGC=  0.033 
VSC=  0.558 
VRC=  0.031 
VKC=  0.018 
VBLD=  0.06 


$' Fraction  liver  tissue 
$' Fraction  fat  tissue,  (0.04  MALE) 

$' Fraction  gut  tissue 
$' Fraction  slow  tissue 
$' Fraction  rapid  tissue 
$' Fraction  kidney  tissue 

$'Est.  fraction  venous  +  arterial  blood  (1+1) 


'PCs  FOR  CH 
CONSTANT  PCHL 
CONSTANT  PCHF 
CONSTANT  PCHS 


PARTITION  COEFFICIENTS 

1.47  $' Liver/blood  partition  coefficient,  CH 

0.48  $'Fat/blood  partition  coefficient,  CH 

1.35  $' Slowly  perfused  tissue/blood  partition,  CH 
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CONSTANT  PCHR  =1.47  $' Richly  perfused  tissue/blood  partition,  CH 

CONSTANT  PCHB  =500.  $'Blood/air  partition  coefficient,  CH 
CONSTANT  PCHG  =  1.47  $' Gut/blood  partition,  CH 

'PCs  FOR  TCOH 

CONSTANT  POHL  =  1.06  $ ' Liver/blood  partition  coefficient,  TCOH 

CONSTANT  POHF  =  1.525  $'Fat/blood  partition  coefficient,  TCOH 
CONSTANT  POHS  =  1.11  .  $' Slowly  perfused  tissue /blood  partition,  TCOH 
CONSTANT  POHR  =  1.06  $' Richly  perfused  tissue/blood  partition,  TCOH 
CONSTANT  POHB  =5000.  $'Blood/air  partition  coefficient,  TCOH 
CONSTANT  POHG  =1.06  $'Gut/blood  partition,  TCOH 
CONSTANT  CLGLUC=0 . 65  $' Renal  clarence  of  GLUC  (L/hr*kg-l) 

CONSTANT  PCBIL  =  70.  $' Bile/blood  partition,  TCOGLUC 

CONSTANT  PCGBO  =  0.31  $' Body /blood  partition,  TCOGLUC 


MOLECULAR  WEIGHTS 


CONSTANT  MWCH  =147.5  $' Molecular  weight  CH  (g/mol) 

CONSTANT  MWTCOH  =149.5  $' Molecular  weight  TCOH  (g/mol) 

CONSTANT  MWGLUC  =325.4  $' Molecular  weight  TCOG  (g/mol) 

CONSTANT  MWTCE  =131.5  $' Molecular  weight  TCE  (g/mol) 

CONSTANT  MWTCA  =163.5  $' Molecular  weight  TCA  (g/mol) 

CONSTANT  MWDCA  =129.  $' Molecular  weight  DCA  (g/mol) 


/ 


/ 


'  EXPOSURE  PARAMETERS  ' 

CONSTANT  CC  =. false .$' Default  to  open  chamber  ' 

CONSTANT  NRATS=  4.  $'Number  of  mice  (for  closed  chamber)  ' 

CONSTANT  KLC  =  0.05  $' First  order  loss  rate  from  closed  chamber  (/hr)' 

CONSTANT  VCHC  =0.75  $' Volume  of  closed  chamber  (L) 

CONSTANT  SODA  =  0.005  $' Volume  of  soda  lime  (L) 

'CH  DOSE 

CONSTANT  PDOSCH=  0.  $'Oral  dose  CH  (mg/kg) 

CONSTANT  KACH  =  1.  $ ' Oral  uptake  rate  CH  (/hr) 

CONSTANT  IVDCH  =0.  $'IV  dose  CH  (mg/kg)  ' 

CONSTANT  COCH  =  0.  $' Concentration  of  inhaled  CH  in  air  (mg/1) 

'TCOH  DOSE 

CONSTANT  PDOSOH= . 0 .  $'Oral  dose  TCOH  (mg/kg) 

CONSTANT  KAOH  =  1.  $ '  Oral  uptake  rate  TCOH  (/hr) 

CONSTANT  KAOHBI=  2.  $ ' Gut  uptake  rate  from  bile  (/hr) 

CONSTANT  IVDOH  =0.  $ ' IV  dose  TCOH  (mg/kg) 

CONSTANT  DOSEG  =0.  $ ' IV  dose  of  glucuronide  (mg) 

' ***  Exposure  definition  ' 

IF  (CC)  RATS  =  NRATS  $' Closed  chamber  simulation  ' 

IF  (CC)  KL  =  KLC 

IF  ( . NOT . CC )  RATS  =  0.  $'Open  chamber  simulation  ' 

IF  (  .  NOT .  CC )  KL  =  0. 

IF  ’(.NOT.CC)  SODA  =  0. 

' (Turn  off  chamber  losses  so  concentration  remains  constant)  ' 

VCH  =  VCHC- RATS *BW- SODA  $'Net  chamber  volume  (L) 

AICH0  =  COCH*VCH*MWCH/24450.  $' Initial  amount  CH  in  chamber  (mg)  ' 

IF  (PDOSCH.EQ.O.)  KACH  =  0.  $' Parenteral  CH  dosing  ' 

IF  (PDOSOH.EQ.O.)  KAOH  =  0.  $' Parenteral  TCO  dosing' 

'  METABOLISM  PARAMETERS 

'CH  METABOLISM 

CONSTANT  RAM  =  0.  $ '  Prof orma  rate  of  TCE->CH  metabolism  ' 

CONSTANT  PCCH  =  0.999  $ ' TCE  that  is  converted  to  CH  (%)E-2 

CONSTANT  PCCOH  =  27.  $'Fst  order  constant  CH->TCOH  (1/hr-lkg) 

CONSTANT  VMCHOC=  0.46  $'Max  Velocity  CH->TCOH  (mg/hr-lkg) 

CONSTANT  KMCHOH=  247.  $ ' Michaelis-Menten  const,  for  CH->TCOH  (mg/L) 
CONSTANT  PCCA  =  0.02  $' Fst . order  constant  CH->TCA  (1/hr-lkg) 

CONSTANT  PCIDA  =  0.3  $' Fst. order  constant  CH->DCA  (1/hr-lkg) 

CONSTANT  KCHC  =2.0  $' first  order  CH  loss-metabol .  by  blood  (/hr/kg)  ' 
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' TCOH  METABOLISM 
CONSTANT  VMTCOC=45 . 0 

CONSTANT  KMTCOH=30 . 

CONSTANT  PCOA  =3.5 
CONSTANT  KOHDC  =  5. 


$'Max.  velocity  of  TCOH->TCOG  (mg/h/kg) 

'DOG  VmTCOC=24 . 6 

$'Michaelis-Menten  const,  for  TCOH->TCOG  (tng/L) 
'DOG  KmTCOH=10 . 

$'Fst. order  constant  TCOH->TCA  (1/hr-lkg) 

$' first  order  OH  loss-metabol .  by  blood  (/hr/kg) 


'***  Timing  commands 

CONSTANT  DAYS  =  2.  $' SIMULATION  DURATION  (days) 

CONSTANT  TCHNG  =  4.0  $' Length  of  inhalation  exposure  (hrs) 

CONSTANT  TINF  =0.01  $'Length  of  IV  infusion  (hrs) 

CONSTANT  POINTS  =  500.  $' Number  of  points  in  plot 

CONSTANT  TGAV  =  0.01  $' Length  of  gavage  infusion  (hr) 

'constant  cint=0.05' 


'***  Scaled  parameters 

QC  =  QCC*BW**0 . 74 
QP  =  QPC*BW**0 . 74 
QL  =  QLC*QC 
QF  =  QFC*QC 
QS  =  0 . 24*QC-QF 
QG  =  QGC*QC 
QR  =  0 . 76*QC-QL-QG 
QK  =  QR-0 . 252*QC 

QB  =  QC- (QL+QG) 
QLB=  QL+QG 
QU  =  QUC*BW 
QBIL=  QBILC*BW 


VL  = 
VF  = 
VS  = 
VG  = 
VR  = 
VK  = 
VB  = 
VBLOOD= 


VLC*BW 
VFC*BW 
0 . 82*BW-VF 
VGC*BW 

0 . 101*BW-VL-VG 
VKC*BW 
BW-VL 
VBLD*BW 


CLGLU= 
VMTCOH  = 
PCOTA= 
PCTCOH= 
VMCHOH= 
PCTCA= 
PIDA= 
IVCHR= 
DOSCH= 
IVOHR= 
DOSOH= 


CLGLUC*BW 
VMTCOC*BW*  *0.7 
PCOA/BW**0 . 3 
PCCOH/BW**0 . 3 
VMCHOC*BW*  *  0 . 7 
PCCA/BW**0 . 3 
PCIDA/BW**0 . 3 
IVDCH*BW/TINF 
PDOSCH*BW 
IVDOH*BW/TINF 
PDOSOH*BW 


KCHD  =  KCHC/BW**0 . 3 
KOHD  =  KOHDC/BW**0 . 3 


$' Values  reset  as  %  of  slow  tissues 

$' Values  reset  as  %  of  rapid  tissues 
$' Fractional  blood  flow  to  kidneys 
'MUST  be  subtracted  from  QR 
$' Blood  flow  through  body  except  liver 
$' Total  flow  of  blood  through  liver (L/h) 
$'Flow  of  urine  (L/hr) 

$ ' Flow  of  bile  (L/hr) 


$' Volume  fat  tissue  (kg)  ' 

$' Values  reset  as  %  of  slow  tissues  ' 

$' Values  reset  as  %  of  rapid  tissues  ' 

$' Volume  of  body  except  liver  (kg)  ' 

$'Est.  volume  of  venous  +  arterial  blood' 

$' Renal  clarence  of  TCOG  (L/hr)  ' 

$' Velocity  TCOH->TCOG  (mg/hr) 

$'Fst. order  TCOH->TCA  (l/h/mouse)  ' 

$'Fst. order  CH->TCOH  (l/h/mouse)  ' 

$' Velocity  2nd  ord.  CH->TCOH  (mg/hr)  ' 
$'Fst. order  CH->TCA  (l/h/mouse)  ' 

$'Fst. order  CH->DCA  (l/h/mouse)  ' 

$' Speed  IV  CH  infusion  (mg/h)  ' 

$'p.o.  dose  CH  per  mouse  (mg)  ' 

$' Speed  IV  TCOH  infusion  (mg/h)  ' 

$'p.o.  dose  TCOH  per  mouse  (mg)  ' 

$ ' Fst  order  loss  rates  (l/h)  per  mouse  ' 


'  repeated  gavage  dosing' 

INTEGER  DAY 

tstop=  days*24. 

CINT  =  tstop/points 

DAY=3 .  $ ' TO  START  GAVAGE  ON  MONDAY  -1,  TUES  0,  WEDN  1,  ETC. 
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END  $'End  of  initial' 

DYNAMIC 


) 

i 


'***  REPEATED  GAVAGE  DOSING 

'GAV  =  FEED  MICE  p.o.  YES=1,  N0=0. 

DISCRETE  CAT1 

INTERVAL  CAT  =  24.  $' EXECUTE  CAT1  EVERY  24  hr 

DAY=DAY+1 

IF (MOD (DAY, 7 ) . GE . 5 )  GOTO  OUT 

GAV  =  1.  $' GAVAGE  =  YES 

SCHEDULE  CAT2  .AT.  T+TGAV  $' SCHEDULE  END  OF  GAVAGE 
OUT . .  CONTINUE 
END  $ ' END  OF  CAT1' 

DISCRETE  CAT2 

GAV  =  0.  $' GAVAGE  =  NO 

END  $ ' END  OF  CAT 2 


ALGORITHM  IALG  =  2  $ ' Gear  method  for  stiff  systems  ' 
'  If  program  hangs -up  at  long  T  with  low  mass  input' 
'  change  to  IALG  =  9  at  .CSL  file  (it  will  use  a  ' 
'  plenty  of  computer  time  to  execute) .  ' 
'  It  may  also  help  to  set  during  the  run  time  -  ' 
'  S  DPSITG= . TRUE .  -  at  . CMD  file 


DERIVATIVE 


/ 


/ 


'  1)  OOPBPK  SUB-MODEL  FOR  CHLORAL  (CH) 

f _ _ _ _ _ _ _ _ _ _ t 


'CICH  =  Concentration  of  CH  in  inhaled  air  (mg/1)  ' 

Cl ZONE  =  RSW ( (T. LT . TCHNG) . OR. CC, 1 . , 0 . ) 

RAICH=  RATS*QP* (CCHA/PCHB-CICH)  -  (KL*AICH) 

AICH=  INTEG (RAICH, AICH0) 

CICH=  AICH/VCH*CIZONE 
CPCH=  CICH*24450./MWCH 
RINHCH=  CICH*QP 
AINHCH=  INTEG (RINHCH, 0 . ) 

' CCHA  =  Concentration  of  CH  in  arterial  blood  (mg/1)  ' 

CCHA  =  (QC*CCHV+QP*CICH) / (QC+ (QP/PCHB) ) 

AUCCHB  =  INTEG ( CCHA , 0. ) 

'ACHX  =  Amount  of  CH  exhaled  (mg)  ' 

CCHX  =  CCHA/PCHB 

CHXPPM  =  (0.7*CCHX+0.3*CICH) *24450. /MWCH 
RACHX  =  QP*CCHX 
ACHX  =  INTEG ( RACHX, 0. ) 

' ACHG  =  AMOUNT  OF  CH  IN  GUT/MOUSE  (mg) 

RACHG  =  QG* (CCHA-CVCHG)  $' PARENTERAL  ADMINISTRATION  ' 

ACHG=  INTEG (RACHG , 0 . ) 

CVCHG=  ACHG/ (VG*PCHG) 

CCHG=  ACHG/VG 

'single  gavage  dosing  ' 

RMRCH  =  -KACH*MRCH 

MRCH  =  DOSCH*EXP(-KACH*T)  $' AMOUNT  REMAINING  IN  STOMACH  (mg) 
RAOCH  =  KACH*MRCH 

AOCH  =  DOS  CH - MRCH  $' TOTAL  MASS  INPUT  FROM  STOMACH  (mg)  ' 

'ACHS  =  Amount  of  CH  in  slowly  perfused  tissues  (mg)  ' 

RACHS  =  QS* (CCHA-CVCHS) 
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ACHS  =  INTEG (RACHS ,  0 .  ) 

CVCHS  =  ACHS/ (VS*PCHS) 

CCHS  =  ACHS /VS 

' ACHR  =  Amount  of  CH  in  rapidly  perfused  tissues  (mg)  ' 

RACHR  =  QR* (CCHA-CVCHR) 

ACHR  =  INTEG (RACHR, 0 . ) 

CVCHR  =  ACHR/ (VR*PCHR) 

CCHR  =  ACHR/VR 

' ACHF  =  Amount  of  CH  in  fat  tissue  (mg)  ' 

RACHF  =  QF* (CCHA-CVCHF) 

ACHF  =  INTEG ( RACHF, 0. ) 

CVCHF  =  ACHF/ (VF*PCHF) 

CCHF  =  ACHF/VF 

'ACHL  =  amount  of  CH  produced  in  liver  per  mouse  ' 

'metabolism  of  TCE-- INTERMEDIATE -->CH.  Production  of  intermediate  ' 

'assumed  to  be  Michaelis-Menten  type.  Decay  of  intermediate  to  yield  ' 
'CH  in  liver  assumed  to  be  a  first  order  reaction.  ' 

RCH  =  ram*PCCH* (MWCH/MWTCE)  $ ' RATE  CONVERSION  TCE->CH  ' 

'ACHL  =  AMOUNT  OF  CH  REMAINING  IN  LIVER/MOUSE  (mg) 

RCHL=QL* (CCHA-CVCHL) +QG* (CVCHG-CVCHL) +RAOCH+RCH-RKCH-RAMCH-RTA-RIDA 
ACHL  =  integ (RCHL, 0 . )  $' AMOUNT  OF  CH  IN  LIVER/MOUSE (mg)  ' 

CVCHL  =  ACHL/ (VL*PCHL)  $' LIVER  VENOUS  BLOOD  CONCENTRATION  CH  ' 
CCHL  =  ACHL/VL  $' LIVER  CONCENTRATION  CH  (mg/L)  ' 

RKCH  =  KCHD*CVCHL*VL  $ ' FST  ORDER  LOSS  RATE-NONSPECIFIC  BIND' 

'RATE  METABOLIZED  TO  TCOH  (mg/hr) :  Michael is -Menten  +  First  order  ' 
RAMCH  = (VMCHOH*CVCHL) / (KMCHOH+CVCHL)  +  PCTCOH* CVCHL* VL 
RTA  =  CVCHL* VL*PCTCA  $ ' RATE  METABOLIZED  TO  TCA  (mg/hr) 

RIDA  =  CVCHL*PIDA*VL  $ ' RATE  METABOLIZED  TO  DCA  (mg/hr) 

AKCH  =  INTEG (RKCH, 0. ) 

AMCH  =  INTEG ( RAMCH, 0. ) 

ATA  =  INTEG (RTA, 0 . ) 

AIDA  =  INTEG (RIDA, 0 . ) 

TOTCH  =  integ (RCH, 0 . )  $' TOTAL  AMOUNT  OF  CH  PRODUCED/MOUSE (mg) ' 

BWCH  =  TOTCH/BW 

MCH  =  BWCH/MWCH  $'MILIMOLES  CH  PRODUCED  (mmoles/kg)  ' 

' IVCH  =  CH  Intravenous  infusion  rate  (mg/hr)  ' 

IVCH  =  IVCHR* (1 . -step (tinf ) ) 

' CCHV  =  Mixed  venous  blood  concentration  CH  (mg/L)  ' 

CCHV  =  (QF*CVCHF+  (QG+QL) *CVCHL+  QS*CVCHS+  QR*CVCHR+  IVCH) /QC 

'BALANCE  OF  CH  IN  MOUSE 

CHBAL = ACHF + ACHL + ACHS +ACHR+ACHG+AMCH+ACHX+ATA+AKCH+AIDA 
'  AMT  IN  MOUSE  (mg) ' 

BLCH  = ( ACHF + ACHL + ACHS  + ACHR +ACHG+ AKCH ) /BW 

' BODY  LEVEL  OF  REMAINING  CH  (MG/KG) ' 
CHIN  =IVCHR*TINF+AINHCH+AOCH+TOTCH  $' DOSE  RECEIVED  (mg/mouse)' 

/ _ > 

'  ***  END  OF  CH  PROGRAM  ***  ' 

/ _ ' 

'  2)  OOPBPK  SUB-MODEL  FOR  TRICHLOROETHANOL  METABOLITE  (TCOH)  ' 

< _ < 

'COHA  =  Concentration  of  TCOH  in  arterial  blood  (mg/1)  ' 

COHA  =  (QC*COHV) / (QC+ (QP/POHB) ) 

AUCOHB  =  INTEG ( COHA, 0. ) 
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'AOHX  =  Amount  of  TCOH  exhaled  (mg) 

COHX  =  COHA/POHB 

OHXPPM  =  (0.7*COHX) *24450. /MWTCOH 
RAOHX  =  QP*COHX 
AOHX  =  INTEG( RAOHX, 0.) 

'AOHG  =  AMOUNT  OF  TCOH  IN  GUT /MOUSE  (mg) 

RAOHG  =  QG* (COHA-CVOHG)  $' PARENTERAL  ADMINISTRATION 

AOHG=  INTEG( RAOHG, 0.) 

CVOHG=  AOHG/ (VG*POHG) 

COHG=  AOHG/VG 

'single  gavage  dosing 
'  MROH  =  AMOUNT  REMAINING  IN  GUTS  (mg) 

MROH  =  DOSOH  *EXP ( -KAOH*T) 

RMROH  =  -KAOH*MROH 
RAOOH  =  KAOH*MROH 

'TOTAL  MASS  INPUT  FROM  STOMACH  (mg) 

AOOH  =  DOSOH  -  MROH 

'  AOHS  =  Amount  of  TCOH  in  slowly  perfused  tissues  (mg) 

RAOHS  =  QS* (COHA-CVOHS) 

AOHS  =  INTEG ( RAOHS, 0. ) 

CVOHS  =  AOHS/ (VS*POHS) 

COHS  =  AOHS /VS 

' AOHR  =  Amount  of  TCOH  in  rapidly  perfused  tissues  (mg) 

RAOHR  =  QR* (COHA-CVOHR) 

AOHR  =  INTEG (RAOHR, 0 . ) 

CVOHR  =  AOHR/ (VR*POHR) 

COHR  =  AOHR/VR 

' AOHF  =  Amount  of  TCOH  in  fat  tissue  (mg) 

RAOHF  =  QF* (COHA-CVOHF) 

AOHF  =  INTEG ( RAOHF, 0.) 

CVOHF  =  AOHF/ (VF*POHF) 

COHF  =  AOHF/VF 

'TOTCOH  =  amt.  of  TCOH  produced  in  liver 
metabolism  of  CH-->TCOH->TCOG.  TCOH  production  in  liver 
assumed  to  be  a  first  order,  and  glucuronidation 
is  a  Michaelis-Menten  type,  not-limited  by  UDPGA  concentration 

RTCOH  =CVCHL*VL*PCTCOH* (MWTCOH/MWCH)  $ ' RATE  CONVERSION  CH->TCOH 

'Rate  TCOH  remaining  in  liver  (mg/hr) 

RTOHL=QL* (COHA-CVOHL) +QG* (CVOHG-CVOHL) +RAOOH+RTCOH-RKTCOH-RGLUC . . . 

+RAGLUC - RAMOH 

ATCOHL=  integ (RTOHL , 0 . )  $' AMOUNT  OF  TCOH  IN  LIVER/MOUSE (mg) 

CVOHL  =  ATCOHL/ (VL*POHL)  $' LIVER  VENOUS  BLOOD  CONCENTRATION  OH 

CTCOHL=  ATCOHL/VL  $' LIVER  CONCENTRATION  TCOH  (mg/L) 

RKTCOH  =  KOHD* CVOHL* VL  $'F-ST  ORDER  LOSS  RATE-NONSPECIFIC  BIND 

RAMOH  =  CVOHL*PCOTA*VL  $ ' RATE  METABOLIZED  TO  TCA  (mg/hr) 

AMOH  =  INTEG (RAMOH, 0 . ) 

AKTCOH=  INTEG (RKTCOH, 0 . ) 


3)  OOPBPK  SUB-MODEL  FOR  GLUCURONIDE  (TCOG)  IN  LIVER,  BILE  AND  BLOOD 


RGLUC  = (VMTCOH* CVOHL) / (KMTCOH+CVOHL)  $' RATE  OF  GLUCURONIDATION 
AGLUC  =  INTEG  ( RGLUC,  0. )  $'AMT.OF  TCOH  GLUCURONIDATED  (mg) 

RG  =  RGLUC*MWGLUC/MWTCOH  $ ' RATE  OF  TCOG  FORMATION  (mg/hr) 
AGA  =  AGLUC*MWGLUC /MWTCOH  $' AMOUNT  OF  TCOG  FORMED  (mg) 


RGL  =  QLB* (CGB-CLVG) +RG-RGBIL 

' AMGL  =  AMOUNT  OF  TCOG  REMAINING  IN  LIVER  (mg) 

AMGL=  INTEG (RGL, 0 . ) +DOSEG  $'AMT.OF  TCOG  REMAIN. IN  LIVER (mg)' 
CLVG=  AMGL/VL  $ ' Pc  BLOOD/TISSUE  =1. 

'Uniform  fast  equilibrium  with  bile  ' 

VBIL  =  QBIL* (T+l) 

CGBIL  =  PCBIL*CLVG 

RGBIL  =  QBIL*CGBIL  $ ' RATE  TCOG  EXCRETION  IN  BILE  (mg/hr)  ' 

AMBIL  =  INTEG (RGBIL , 0. )  $' AMOUNT  OF  TCOG  IN  BILE  (mg) 

ATGB  =  AMBIL*MWTCOH/MWGLUC  $' AMOUNT  OF  TCOH  FROM  TCOG  IN  BILE  ' 

'  MRGLUC=  AMOUNT  OF  TCOH  FROM  HYDROLYZED  TCOG  REMAINING  IN  GUTS  (mg)  ' 

MRGLUC=  ATGB*EXP(-KAOHBI*T) 

RAGLUC=  KAOHB I *MRGLUC  $ ' RATE  OF  ABSORPTION  OF  HYDROLYZED  TCOG' 
AGABS  =  ATGB  -  MRGLUC  $ ' TCOG  INPUT  FROM  GUTS  (mg) 

TCOHM  =  AMOH  +  AGLUC  -  AGABS  $' TOTAL  AMOUNT  OF  TCOH  METABOLIZED' 


TOTCOH=  integ (RTCOH , 0 . ) 
bwtcOH=  totCOH/bw 

MOTCOH=  BWTCOH/MWTCOH  $'MILIMOLES  TCOH  PRODUCED  (mmoles/kg)  ' 


'Fast  uniform  diffusion  within  the  volume  of  distribution  ' 

RBODYG=  QB* (CGB-CBVG) -RGU  $ ' RATE  OF  TCOG  REMAINING  IN  BODY (mg/h) ' 
ABODYG=  INTEG ( RBODYG , 0. ) 

CBOG=  ABODYG/VB  $' TISSUE  CONCENTRATION  OF  TCOG  ' 

CBVG=  CBOG/PCGBO  $'CONC.TCOG  IN  ULTRAFILTR . BLOOD  ' 

CGB=  (QB*CBVG  +  QLB*CLVG) /QC  $' CONCENTRATION  OF  TCOG  IN  BLOOD' 
AGB=  INTEG ( CGB,  0. )  $' AMOUNT  OF  TCOG  IN  BLOOD  (mg) 

I _ _  _  __  _ _ _ _ _  _ _ _ _ _ _ _ _ . _ , _ _ _ / 


'  AGU  =  AMOUNT  OF  TCOG  IN  URINE  (mg)  ' 

VU  =  QU* (T+l) 

RGU  =  CBVG*  CLGLU  $ '  RATE  OF  TCOG  REMOVAL  IN  URINE  (mg/hr)' 

AGU  =  INTEG (RGU, 0 . )  $' AMOUNT  OF  TCOG  EXCRETED  IN  URINE  (mg)' 

CGUEX  =  AGU/VU  $ ' CUMULATED  CONC .  TCOG  EXCRETED  IN  URINE  (mg/L) ' 
CGU  =  RGU/QU  $' CONCENTRATION  OF  TCOG  IN  URINE  (mg/L)' 


GLUBAL=  AGA  +  DOSEG  -  (AMGL  +  AMBIL  +  ABODYG  +  AGU) 


***  CONTINUE  WITH  COMMON  PART  FOR  TCOH  AND  TCOG  *** 


' IVOH  =  TCOH  Intravenous  infusion  rate  (mg/hr) 

IVOH  =  IVOHR* (1. -step (tinf) ) 

' COHV  =  Mixed  venous  blood  concentration  TCOH  (mg/L) 

COHV  =  (QF*CVOHF+  (QG+QL) *CVOHL+  QS*CVOHS+  QR*CVOHR+  IVOH) /QC 


'BALANCE  OF  TCOH  IN  MOUSE 


'  AMOUNT  OF  TCOH  IN  MOUSE  (mg) ' 

OHBAL  =  AOHF+ATCOHL+AOHS+AOHR+AOHG+AMOH+AOHX+AGLUC+AKTCOH- AGABS 
BLOH  =  ( AOHF + ATCOHL + AOHR + AOHS  + AOHG+ AKTCOH+AGLUC - AGAB S ) / BW 

BODY  LEVEL  OF  REMAINING  AND  EXCRETED  TCOH  (MG/KG) ' 
BLOHF  =  BLOH+ ( AGABS /BW) - (AGLUC/BW) 

BODY  LEVEL  OF  FREE  -  NON  GLUC .  TCOH  (MG/KG)' 
OHIN  =  IVOHR*TINF+AOOH+TOTCOH  $ ' DOSE  RECEIVED/MOUSE  (mg) 


***  END  OF  TCOH  AND  TCOG  PROGRAMS  *** 


*  RATES  AND  AMOUNTS  OF  ACID  METABOLITES  * 


RTCA  =CVCHL*VL*PCTCA* (MWTCA/MWCH)  $ 'RATE  CONVERSION  CH->TCA  ' 


'Additional  amount  of  TCA  is  generated  from  TCOH  (rate  =  RAMOH) : 
RCAOH  =  RAMOH* (MWTCA/MWTCOH) 


RCHDC=RIDA* (MWDCA/MWCH)  $ ' RATE  CONVERSION  OF  CH->DCA  (mg/h) ' 


ATCA  =  INTEG  ((RTCA  +  RCAOH),  0.)  $' Amount  of  TCA  (mg)' 

ADCA  =  INTEG  (RCHDC,  0.)  $' Amount  of  DCA  (mg)' 


TERMT  (T.GE.TSTOP) 

END  $'End  of  derivative 

END  $'End  of  dynamic 

END  $ ' End  of  program 


'***************************************************************'*******' 


*.CMD  File 


Cabana  and  Gessner  (1970)  Mice: 


PROCED  Fig4a 

prepar  t , BLCH , cchv , cchl 

'male  mouse' 

SET  TITLE= ' CH  iv  500mg/kg  MOUSE (Cabana, 1970) ' 
s  symcpl= . true . ,  ivdch=500 . , PCCOH=15 . , kohdc=7 . 9 
s  tchng=0, days=0 . 25 , qlc= . 24, VMCHOC=3 . 23 , kmtcoh=l . 
s  qpc=30, qcc=16 . 5 , vfc= . 04 , qf  c= . 05 , kchc=l . 04 , vmtcoc=15 
s  bw= .  03  ,grdcpl= .  f . ,  wesitg= .  f .  KMCHOH=0 .0675 

DATA 

T  BLCH 

0.  500.0  INITIAL 

0.08  475.0 

0.2  280.0 

0.25  230.0 

0.333  200.0 

0.5  150.0 

0.75  90.0 

1.  60.0 

1.5  48.0 

2.  34.0 

3.  50.0 

4.  44.0 

6.  48.0 

END 

START 

PLOT  BLCH, ' log' , ' lo' =10 . , 'xhi' =6 . 

END 

PROCED  Fig4b 

prepar  t,BLOH, cchv, cchl, cohv 
'male  mouse' 

SET  TITLE= ' CH  iv  500mg/kg  MOUSE  (Cabana, 1970) ' 
s  symcpl= . true . ,  ivdch=500 . ,pccoh=15 , kohdc=7 . 9 
s  tchng=0, days=0 . 25 , qlc= . 24 , vmchoc=3 .23 , kmtcoh=l. 
s  qpc=30 , qcc=16 . 5 , vfc= . 04 , qfc= . 05, kmchoh=0 . 0675, 
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s  bw= . 03 , grdcpl= .  f . , wesitg= .  f . , kchc=l . 04 , vmtcoc=15 . 

DATA 

T  BLOH 

0.  0.00001  INITIAL 

0.07  65.0 

0.15  120.0 

0.25  135.0  > 

0.33  170.0 

0.5  200.0  * 

0.75  230.0 

1.  235.0 

1.5  250.0 

2.  255.0 

3.  250.0 

4.  260.0 

6.  252.0 

END 

START 

PLOT  BLOH, 'log' , ' lo'=10. , 'xhi'=6. 

END 

PROCED  Fig4c 

prepar  t, BLOHF, cchv, cchl , cohv 
'male  mouse' 

SET  TITLE= ' CH  iv  500mg/kg  MOUSE (Cabana, 1970) ' 
s  symcpl= . true . ,  ivdch=500 . ,pccoh=15 . , kmtcoh=l . 
s  tchng=0 , days=0 . 25, qlc= . 24 , vmchoc=3 . 23 ,kohdc=7 . 9 
s  qpc=30 , qcc=16 . 5 , vf c= . 04 , qfc= . 05 , kmchoh=0 . 0675, 
s  bw= . 03 , grdcpl= . f . , wesitg= . f . , kchc=l . 04 , vmtcoc=15 . 

DATA 

T  BLOHF 

0.  0.00001  INITIAL 

0.063  60.0 

0.125  110.0 

0.25  130.0 

0.3  155.0 

0.45  190.0 

0.7  225.0 

1.  225.0 

1.5  215.0 

2.  190.0 

3.  150.0 

4.  105.0 

6.  90.0 

END 

START 

PLOT  BLOHF, 'log' , 'lo'=10. , 'xhi'=6. 

END 

- - END  OF  MICE  DATA - ' 

' - Garrett  and  Lambert  (1993)  Dogs: - ' 

PROCED  Fig 5 a 
prepar  t , cohv, cgb 

SET  TITLE= ' TCOH  iv  30mg/kg  DOG  (Garrett, 1973) ' 
s  bw=2  0 . ,  ivdoh=30 . , days=0 . 4 , qcc=18 . , qpc=20 . , kchc=l .  04 
s  ivdch=0 . , vmtcoc=15 ,  kmtcoh=l . , pccoh=15 , kmchoh=0 . 0675 
s  symcpl= . true . , grdcpl= . f . , wesitg= . f . , vmchoc=3 .23 
s  kohdc=7 . 9 , doseg=0 . 

% 

DATA 
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T 

COHV 

0.0001 

100.8 

0.05 

67.2 

0.26 

34.9 

0.42 

21.8 

0.94 

15.1 

1.3 

10.9 

2.24 

6.1 

2.76 

4.9 

3.39 

3.3 

4.17 

1.8 

4.7 

1.4 

5.8 

1.1 

7.7 

END 

START 

0.7 

PLOT  COHV, 
END 

' lo' =0 

'xhi' =8 . 


PROCED  Fig5b 
prepar  t , cohv , cgb 

SET  TITLE= ' TCOH  iv  60mg/kg  DOG  (Garrett, 1973) ' 
s  bw=10.,  ivdoh=60 . , days=0 . 2 , gcc=18 . , qpc=20 . , kchc=l . 04 
s  vmtcoc=15,  kmchoh=0.0675,kmtcoh=l. ,vmchoc=3 .23,doseg=0. 
s  symcpl= . true . ,grdcpl=.f . ,wesitg=.f . ,pccoh=15 , kohdc=7 . 9 


DATA 

T 

COHV 

0.0001 

109.2 

0.05 

102.1 

0.1 

97.4 

0.15 

88.2 

0.2 

73.9 

0.3 

53.1 

0.5 

30.2 

0.75 

21.4 

1. 

16. 

1.26 

14.3 

1.5 

11. 

1.75 

12.6 

2  . 

6.7 

2.5 

4.7 

3  . 

3.6 

3.5 

2.7 

4. 

1.5 

END 

START 

PLOT  COHV, 
END 

' hi ' =120 

=4 . 


PROCED  FigSa 
prepar  t , cohv , cgb 

SET  TITLE= ' TCOG  iv  450mg/kg  DOG (Garrett, 1973) ' 
s  bw=ll . 85, doseg=450 . , days=0 . 21 , ivdoh=0 . ,pcgbo=0 . 311 
s  gcc=18 . , gpc=20 . , vmtcoc=15 , kmtcoh=l . , pccoh=15 , kchc=l . 04 
s  symcpl= . true . , wesitg= . f . ,grdcpl= . f .  ,  vmchoc=3 . 23 ,kmchoh=0 .00675 

DATA 

T  CGB 

0.0001  183. 

0.02  146.4 

0.07  119. 
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0.117 

91.5 

0.14 

82.4 

0.16 

73.2 

0.28 

54.9 

0.33 

43.9 

0.42 

36.6 

0.51 

31.1 

0.7 

23 . 8 

0.89 

15.9 

1.05 

11.7 

1.21 

7.7 

1.38 

6.4 

1.56 

4.2 

1.8 

3.3 

2.1 

2.3 

2.4 

1.65 

2.7 

0.95 

3. 

0.88 

3.5 

0.49 

4. 

0.44 

4.5 

END 

START 

0.23 

PLOT  CGB, 
END 

' xhi ' =4 . 

PROCED  Fig6b 

prepar  t, cohv, cgb, cgbil 

SET  TITLE= ' TCOH  iv  llmg/kg  DOG (Garrett,  1973)  ' 
s  bw=13 . ,  doseg=0.(  ivdoh=ll . , days=0 . 2 ,kchc=l . 04 
s  vmtcoc=15,  kmtcoh=l . , vmchoc=3 . 23 ,pccoh=15 , kohdc=7 . 9 
s  symcpl=.t.,  grdcpl= . f . ,  wesitg= .  f .  ,ktnchoh=0 . 0675 

DATA 


T 

COHV 

CGB 

CGBIL 

0.05 

19.3 

. 

. 

0.1 

15.1 

. 

. 

0.2 

10.1 

5.9 

. 

0.5 

5.9 

13.7 

. 

0.75 

4.6 

15.6 

. 

1. 

3.7 

16.5 

823.5 

1.25 

2.9 

16 . 7 

1.5 

2.3 

14.6 

. 

1.75 

1.9 

11.9 

. 

2  . 

1.85 

10.1 

585 . 6 

2.5 

1. 

6.8 

. 

3. 

0.8 

4.6 

179.3 

3.5 

3.8 

4. 

2.2 

102.5 

END 

START 

PLOT  COHV, 

& 

H- 

II 

to 

O 

,  CGB ,  i 

CGBIL, 'xhi 

END 

' - END  OF  DOGS  DATA . . . 

' - CH  and  TCOH  IV  Mice  OUR  data - - 

PROCED  Fig7a 

prepar  t, cchv, cchl , clvg, cgb 
'male  mouse' 

SET  TITLE= ' CH  iv  lO&lOOmg/kg  MOUSE' 

s  symcpl= .  true . , ivdch=100 . , ivdoh=0 . , POHG=l . 06 , kchc=l . 04 , kohdc=7 . 9 
s  tchng=0 , days  =  0 . 0042 ,  glc= . 24 , PCHS=1 . 226 , KMCHOH=0 . 0675 , kmtcoh=l . 

S  PCHR=1 . 3  98 , PCHG=1 . 47 , POHL=l . 06 , POHF=l . 525 , POHS=l . 11 , POHR=l . 06 , 


38 


s  gpc=3 0 ,  qcc=16 . 5 , vf c= . 04 , qf c= . 05 , tinf =0 . 0125 , PCCOH=15 , VMCHOC=3 .23 
s  bw= .  042 , grdcpl= .  f .  , wesitg= .  f .  , PCHL=1 .398, PCHF=0 . 366 , vmtcoc=15 . 


DATA 

T 

cchv 

ivdch 

0. 

0. 

100. 

0.0149 

304.62 

0.0149 

265.78 

0.0149 

232 . 88 

0.0158 

427.5 

0.0485 

95.28 

0.0485 

69.36 

0.0485 

70.91 

0.0485 

111.87 

0.09 

56.03 

0.09 

30.18 

0.09 

62.17 

0.09 

45.06 

0.097 

52.01 

0. 

0. 

10. 

0.0143 

24.26 

0.0143 

12.60 

0.0143 

8.85 

0.0143 

11.79 

0.0476 

3.26 

0.0476 

4.75 

0.0476 

5.60 

0.0476 

3.30 

0.088 

0.31 

0.088 

0.75 

0.088 

0.40 

END 

START  smooth 
PLOT  cchv,  ' xhi ' 

=  0.1 

END 


INITIAL 


INITIAL 


PROCED  Fig7b 

prepar  t, cohv, BLOH, BLCH, cgb, ctcohl,  clvg 
'male  mouse' 

SET  TITLE= ' TCOH  iv  10&100mg/kg  MOUSE' 

s  symcpl= . true . , ivdch=0 . , ivdoh=100 . , POHG=l . 06 , vmchoc=3 . 23 , 
s  tchng=0, days=0 . 0042 , qlc= . 24 , PCHS  =  1 . 35, pccoh=15, kmchoh=0 . 0675 
S  PCHR=1 . 47 , PCHG=1 . 47 , POHL=l . 06 , POHF=l .525, POHS  =  l . 11 , POHR=l .  06  , 
s  qpc=30 , qcc=16 . 5 , vfc= . 04 , qfc= . 05, tinf=0 . 0125 , vmtcoc=15 . , kchc=l . 04 
s  bw= . 042 , grdcpl= . f . , wesitg= . f . , PCHL=1 . 47 , PCHF=0 . 48 , kmtcoh=l . 
s  kohdc=7.9 


DATA 

T 

COHV 

0. 

0. 

0.012 

359.1 

0.012 

187.1 

0.012 

200 . 7 

0.0197 

198.9 

0.0458 

78.56 

0.04945 

50.43 

0 . 04945 

67 . 06 

0.04945 

52 . 79 

0.0892 

40.28 

0.0892 

48.88 

0.0892 

53.56 

0. 

. 

ivdoh 

100.  INITIAL 
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INITIAL 
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0.015 

14.99 

0.0162 

9.98 

0.0162 

10.39 

0.0162 

12.22 

0 . 0484 

4.47 

0 . 0484 

3.09 

0 . 0484 

3.90 

0.0484 

4.36 

0.0847 

2.70 

0.0860 

2.46 

0.0886 

2 . 59 

0.0886 

4.85 

0.0886 

1.83 

END 

START  smooth 

PLOT  cohv,  ' xhi '=0.1 

END 

' . . - . END  OF  OUR  MICE  DATA 
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